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Abstract: The complete tensor structure of the quark-gluon vertex in Landau gauge is 
determined at two kinematical points ('asymmetric' and 'symmetric') from lattice QCD in 
the quenched approximation. The simulations are carried out at /3 = 6.0, using a mean- 
field improved Sheikholeslami-Wohlert fermion action, with two quark masses ~ 60 and 
115 MeV. We find substantial deviations from the abelian form at the asymmetric point. 
The mass dependence is found to be negligible. At the symmetric point, the form factor 
related to the chromomagnetic moment is determined and found to contribute significantly 
to the infrared interaction strength. 
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1. Introduction 

The quark-gluon vertex describes the coupling between quarks and gluons, and is thus 
one of the fundamental quantities of QCD. In perturbation theory, a complete calculation 
has been performed to one loop [Q, and partial two- and three- loop calculations have 
been performed for specific gauges and kinematics f|. Nonperturbatively, however, it 
remains largely unknown. In [|, || the first steps were taken towards a nonperturbative 
determination, by way of a quenched lattice calculation of the form factor containing the 
running coupling in two different kinematics in the Landau gauge. 

The Dyson-Schwinger equation (DSE) for the quark propagator contains the quark- 
gluon vertex, and normal practice has been to truncate the hierarchy of DSEs by providing 
an ansatz for the vertex. However, if a realistic gluon propagator, obtained from the 
coupled ghost-gluon(-quark) DSEs jj], ||, |9) and consistent with lattice data flO], |ll], |l2| is 
used, dynamical chiral symmetry breaking appears to be quite sensitive to the details of 
the ansdtze employed |J. It therefore appears highly desirable to obtain 'hard' information 
about the full infrared structure, not only the part containing the running coupling. 

In this paper we take the first steps towards this aim, by determining all the nonzero 
form factors at the two kinematic points used in ||, namely q = and q = —2p, where q 
is the gluon momentum and p is the momentum of the outgoing quark leg. At the same 
time we also study the quark mass dependence by using two different quark masses for the 
vertex at q = 0. Some preliminary results have already been presented in [13]. 

The quark-gluon vertex is related to the ghost self-energy through the Slavnov-Taylor 
identity, 

<7 M r>, q) = G{q 2 ) [(1 - B(q,p + q^^p) - S~Hp + ?)(!- B(q,p + q))] , (1.1) 
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where G(q 2 ) is the ghost renormalisation function and B(q, k) is the ghost-quark scattering 
kernel. Evidence from lattice simulations [14] and Dyson-Schwinger equation studies M, M, 
|9| indicate that G(p 2 ) is strongly infrared enhanced, and this should also show up in the 
quark-gluon vertex. On the other hand, nontrivial structure in the ghost-quark scattering 
kernel, which has usually been assumed to be small, may also be realised in the vertex. 

The rest of the paper is structured as follows: In section ^ we briefly present our 
notation and procedure, referring to || for the details. In section || we present results for 
the vertex at the asymmetric point and compare to the abelian (quark-photon) vertex, 
which is completely determined by the Ward-Takahashi identity at this point. In section 
||]we present results for the vertex at the symmetric point, including the 'chromomagnetic' 
form factor tq. Finally, in section || we summarise our results and discuss prospects for 
further work. Some tree-level formulae used in the analysis are given in the Appendix. 



2. Notation and procedure 



Throughout this article, we will be using the same notation as in [0], and we refer to that 
article for a detailed discussion of our notation and procedure. We write the one-particle 
irreducible (proper) vertex (see fig. [[]) as A^(p,q) = t a A^(p,q), where p and q are the 
outgoing quark and gluon momentum respectively. The incoming quark momentum is 
denoted k. 

We will be operating in the Landau gauge, where, 
as discussed in ||, only the transverse-projected part 
of the vertex can be studied away from q = 0. We 
will therefore define the transverse-projected vertex 
as 



mi a 
a 



9*8 



A P (p,q)=P^{q)A v {p,q) 



W ~ ^-r) A u(p,q) ■ 
(2.1) 

In a general kinematics the vertex can be decomposed 
into 12 independent vectors which we can write in 
terms of vectors Li,Tj and scalar functions \i,Ti as 
described in ||: 




+ g 



Figure 1: The quark-gluon vertex. 



q) = -ig k2 ) L iAPi ?) ~ n{p 2 , q 2 , k 2 )Ti^{p, q) . (2.2) 



i=i 



i=i 



We will here be focusing on the two specific kinematics defined in || and related there 
to the MOM and MOM renormalisation schemes — namely, the 'asymmetric' point q = 
(i.e., p 2 = k 2 ;q 2 = 0) and the 'symmetric' point q = —2p (i.e., p 2 = k 2 = q 2 /4). In the 
asymmetric kinematics, the vertex reduces to 



A^O) =-ig Xi(p 2 ,0,p 2 h, 



4\ 2 (p 2 ,0,p 2 ) 



2i\ 3 (p 2 ,0,p 2 )p, 



(2.3) 
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while in the symmetric kinematics we have 



A M (-g/2,g) 



A£(-<?/2,g) 



Ai(g 2 /4, q 2 , q 2 /4h» + Mq 2 /A, q\ q 2 /4) (fa - q 2 ^) 
- ir 5 (q 2 /4,q 2 ,q 2 j4)a liv q v 



W/4,gW/4)( 7M - iq^lq 2 ) - ir 5 (q 2 /4,q 2 ,q 2 /4)a t 



[LVHV 



(2.4) 



(2.5) 



where on the last line, in the transverse-projected vertex, we have written X[ = Ai — q 2 r%. 

In an abelian theory (QED), the Ward-Takahashi identities imply that the form factors 
\{i = 1, 2, 3) are given uniquely in terms of the fermion propagator, 



S(p) 



Zip 2 



i i>A(p 2 ) + B[p 2 ) i i> + M(p 2 ) 
In the kinematics we are considering, they are given by 



A« E V,o, P 2 ) = ~Mp 2 ); a« e V,o,p 2 ) = -^ 2 B{ P 2 ) 



QED, 2 



(2.6) 

(2.7) 
(2.8) 



The deviation of the QCD form factors from these expressions thus give us a measure of 
the purely nonabelian nature of the theory. Note that A4, which is identically zero in QED, 
is zero also in QCD at these two particular kinematic points. 

The bare (unrenormalised) quantities \3,t§ and \± (at the symmetric point) can be 
obtained by tracing the lattice with the appropriate Dirac matrix (the identity, a^ u and 
7 /t respectively): 



A 3 (p 2 ,0,p 2 ) = ^^p^ReTrA^M); 
r 5 (g 2 /4, q 2 , q 2 /A) = ~ ? „ J- ReTr a^(-q/2, q) ; 



490 



Ai(^ 2 /4,^,g74) 



2 „2 



^E^- Imlv > A ?(-«A?)- 



(2.9) 
(2.10) 
(2.11) 



At the asymmetric point, Ai and A2 both come with the same Dirac structure. To separate 
them, we first determine Ai as described in || by setting the 'longitudinal' momentum 
component p^ to zero, and then obtain A2 by 

A 2 (p 2 ,0,/) = ^^(^-ImTr 7M A>,0) + A^p 2 , 0,p 2 )) . (2.12) 



In order to make the lattice form factors more continuum-like, we employ tree-level 
correction, as discussed in [11, [y|. The tree-level correction of Ai is described in Q, 
although at the symmetric point we have here refined the correction procedure, as described 
in appendix [A|. In the case of A2,A3 and T5, these are simply zero at tree level in the 
continuum, while they are non-zero on the lattice with the action and parameters we are 
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using. We therefore have to subtract off the lattice tree-level forms. The details of this are 
given in appendix |A[ Not unexpectedly, this procedure leads to large cancellations which 
make our results unreliable at large momenta. 

As always, the quantities obtained from the lattice are bare (unrenormalised) quanti- 
ties. The relation between renormalised and bare quantities is given by 

< = Z 2 1/2 V>; / = Z 2 1/ V; Al = zl /2 A„; g = Z g g; £o = Z 3 £ , (2.13) 

where Z2,Z^,Z g are the quark, gluon and vertex (coupling) renormalisation constants 
respectively. The renormalised quark and gluon propagator and quark-gluon vertex are 
related to their bare counterparts according to 

5 bare (p; a) = Z 2 {n\ a)S(p; /i) ; D h ™ c (q 2 ; a) = Z 3 (jjl; a)D(q 2 ; M ) ; (2.14) 
A£ arc 0p, q; a) = Z${p, a)A M (p, q; fj) . (2.15) 

Renormalisation may be carried out in a momentum subtraction scheme. For the quan- 
tities computed at the asymmetric point, we will use the MOM scheme defined in || 
requiring that Ai(/i 2 , 0, fi 2 ) = 1; while for the quantities at the symmetric point we will 
use a modification of the MOM scheme, requiring A / 1 (/i 2 /4, /i 2 , /i 2 /4) = 1. In both cases we 
choose [i = 2 GeV as our renormalisation scale. We can then easily match our results on 
to perturbation theory in the ultraviolet, using the associated (MOM or MOM) running 
coupling. 

We use the same ensemble and parameters as in [p|. The Wilson gauge action is used 
at (3 = 6.0 on a 16 3 x 48 lattice. The Sommer scale provides an inverse lattice spacing of 
2.12 GeV. The mean-field improved SW action is adopted with off-shell improvement in the 
associated propagators. Further details may be found in ||. In order to study the quark 
mass dependence of the vertex, we have used two values for the hopping parameter, k = 
0.137 and 0.1381, corresponding to a bare quark mass m ~ 115 and 60 MeV respectively. 



3. Asymmetric point 



First, we investigate the mass dependence of the Ai form factor, which was already studied 
in p). Since in this paper we are primarily concerned with the deviation from the abelian 
(Ball-Chiu) form, we show, in figure 0, the quantity Z{p 2 )X\{p 2 , 0,p 2 ), which in an abelian 
theory would be constant. The clear infrared enhancement observed in Q is confirmed, and 
we also see that the mass dependence of this quantity is negligible. The slight difference 



in Ai between the two masses observed in [13] is in other words entirely due to the mass- 
dependence of the quark renormalisation function. 

we have fitted the tree- 



In order to compare our results with the abelian forms (2. 



level corrected quark propagator [15] to the following functional forms [16], 

Z(p 2 ) = l/A(p 2 ' ' 



k\ 1 



aM{p 2 ) = aB{p 2 )/A{p 2 



a 2p2 _|_ p 
,2(a-l) 



(aV) a + e 



+ m f 



(3.1) 



(3.2) 
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Figure 2: The unrenormalised form factor Xi(p 2 ,Q,p 2 ) multiplied by the quark renormalisation 
function Z(p 2 ), as a function of p. In an abelian theory, this would be a p-independent constant. 



where k,c,l,c m ,l m ,a and m are fit parameters. The best fit values are given in table |. 
When comparing with the renormalised vertex we use the values obtained from the quark 
propagator renormalised at 2 GeV, which amounts to dividing the unrenormalised values 
by Z(4GeV 2 ). From these fits, we can then derive the abelian form factors (El 



m (GeV) 


k 


c 


I 


Cm 




a 


m f 


60 


1.075 


0.218 


0.326 


0.0261 


0.400 


1.232 


0.0258 


115 


1.045 


0.208 


0.316 


0.0357 


0.484 


1.361 


0.0670 



Table 1: Fit parameters for best fits of the quark propagator to the functional forms (3.1) and 



(3.2). All fits have been performed to data surviving a cylinder cut with radius 1 unit of spatial 
momentum, up to a maximum momentum of pa = 1.2 for the lighter quark mass and 1.4 for the 
heavier mass. 



We will also compare our results with the one-loop Euclidean-space expressions, 
A?V,0,p 2 ) 



f 1 

16vr 2 Ap 2 



m 



P 



+ ^lnfi + 4) MaV + d-OO] 



g 2 m 



16-7T 2 p 2 



(3 + C)C F - (3 + 20 



C A 



m 2 , / p 2 
-s- In 1 + - 




(3.3) 



(3.4) 



where the group factors C F = | and Ca = 3 in QCD, and the gauge parameter £ = 
in Landau gauge. In order to match this to our lattice results, we renormalise both the 
lattice and perturbative data in the MOM scheme. From the data of fig. 6 in we 



find that l/Z^P M (2GeV, a) = 1.391? at to = 115 MeV. From this we determine the 



i+6 
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renormalised form factors A23 



z MOM A iat The M q M i_i oop values are determined by 



using Qrr^rz (2GeV) 
to MOM v ' 



2.21(10) and multiplying by 



evaluating the expressions ( |3,3| ), ( pT4| ) 
Z MOM/ z m = ! 069, obtained from eq. (7.2) of §. 

In figure ^ we show the form factor A2 as a function of p, for the heavier quark mass. 
We see that it is greatly enhanced both compared to the Ball-Chiu form ( |2.8| ) and the one- 
loop form ( |3.3| ), and only approaches these around or above 3 GeV. 1 In figure || we show 
the dimensionless quantity 4p 2 \2(p 2 , 0,p 2 ) as a function of p. This quantity measures the 
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Figure 3: The renormalised form factor A2(p 2 ,0,p 2 ) as a function of p. Also shown is the abelian 



(Ball-Chiu) form of (|2.8[) and the one- loop form of (3.3) 




Figure 4: The renormalised form factor 4p 2 \2(p 2 , 0,p 2 ) as a function of p. 



1 In Il3] there was an error of a factor of 4 in the normalisation of A2, which gave the false impression 
that our numerical results agree almost perfectly with the Ball-Chiu form. 
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relative strength of this component compared to the tree-level Ai. We see that A2 becomes 
comparable in strength to Ai for the most infrared points. 

In figure [| we show \%{p 2 , 0,p 2 ) as a function of p. Here we have performed a 'cylinder 



cut' [17] with radius 1 unit of spatial momentum to select data close to the 4-dimensional 
diagonal. We see that it coincides within errors with the Ball~Chiu form ( |2.8D , and ap- 
proaches the one-loop form at about 2 GeV. We also see that the quark mass dependence 
of both A2 and A3 is very weak. A3 becomes somewhat larger as the quark mass is 



• m=115MeV 

O m = 60 MeV 

BC (115 MeV) 

BC(60MeV) 

— - 1 -loop (115 MeV) 




p (GeV) 



Figure 5: The renormalised form factor As(p 2 ,0,p 2 ) as a function of p. Also shown is the abelian 



(Ball-Chiu) form of (|2.8|) and the one-loop form of (3.4) 
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Figure 6: The renormalised form factor \^{p 2 ,0,p 2 ) multiplied by twice the quark momentum 2p, 
as a function of p, for m = 115 MeV. This dimensionless quantity gives a measure of the relative 
strength of A3. 



-7- 



decreased, which corresponds to the effect of dynamical chiral symmetry breaking being 
relatively larger for a smaller bare mass. 

In figure ^ we show 2p\s(p 2 , 0, p 2 ) as a function of p. This quantity is dimensionless and 
measures the relative strength of A3 compared to the tree- level X\. For the most infrared 
points, A3 can also be seen to contribute significantly to the interaction strength, although 
clearly not as strongly as Ai and A2- 

In order to see more clearly the relative strength of all three components of the vertex, 
in figure [?] we show the dimensionless quantities Ai,4p 2 A2 and 2p\^ for the heavier quark 
mass. In this figure, the hierarchy of strengths Ai > A2 > A3 is evident. 
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Figure 7: The dimensionless form factors Ai, 4p 2 A2 and 2p\3 at the asymmetric point, as a function 
of p, for m = 115 MeV. 



4. Symmetric point 

Since we have already established that the dependence of the vertex on the quark mass 
is very weak, in this section we will only be using one quark mass, m ~ 115 MeV. We 
will also in this section make use of the lattice momentum variables K{p) = sm{pa)/a 
and Q(q) = 2 sin(qa/2)/a = —2K{p). These momentum variables appear in the lattice 
treedevel expressions for the form factors we will be studying, as well as in the transverse 
projector, and are thus appropriate variables to use. 

In figure ^ we show A' x at the symmetric point as a function of |Q(g)|. In contrast 
to in ||, the treedevel correction here has been carried out on each Lorentz component 
of the vertex separately, as explained in the Appendix. These results should therefore 
be more reliable than those shown in ||. We have also performed a cylinder cut on the 
data with a radius of 2 units of spatial momentum in q. From these data, we deter- 
mine l/ZMP M (2GeV,a) = 0.95(8). Multiplying by Z 2 Z^ /2 determined in @ we also find 
5 MUM (2GeV) = 1.47(15). The ratio of renormalisation constants is ZfP M /Z]f = 1.093. 
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Q(q) = 2K(p) (GeV) 

Figure 8: The unrenormalised form factor X[(p 2 , q 2 ,p 2 ) at the symmetric point q = —2p, as a 
function of the gluon momentum q. The data shown are those surviving a cylinder cut with radius 
2 units of spatial momentum in q. 
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Q(q) = 2K(p) (GeV) 

Figure 9: The renormalised form factor X[(p 2 , q 2 ,p 2 ) at the symmetric point q = —2p, as a function 
of the gluon momentum q. Also shown is the one-loop form from B. 



This is used to determine the one- loop A' x , shown together with the renormalised lattice 
in figure |j| 

In figure |l0| we show the form factor t§ as a function of the gluon momentum q. The 
same cylinder cut has been performed as in fig. ||. We see that, although t§ is power 
suppressed in the ultraviolet, it rises very significantly for q < 2 GeV. Although this 
form factor is related to the chromomagnetic moment, and as such is expected to be of 
phenomenological importance, it has not previously been included in QED-inspired model 
vertices commonly used in, e.g., DSE-based studies. However, work is in progress to provide 



% o.i 

a 



0.01 



0.001 



Hi 



J I I L_ 

Q(q) = 2K(p) (GeV) 



10 



Figure 10: The renormalised form factor T5 at the symmetric point as a function of the gluon 
momentum q. The data shown are those surviving a cylinder cut with radius 2 units of spatial 



momentum in q. Also shown is the one-loop form of (4.1) 



an analytical, nonperturbative expression for this and the other form factors in the purely 
transverse part of the vertex [jl8|. We will also compare our lattice results to the one-loop 
T5, which in Euclidean space is given by 



m 



16vr 2 12s 2 



(i-O 



8Cf + £Ca - Ca 



m 2 (l - f) 



- (2C F - C A ) 



+ C a - 



4s 2 (1 







s 2 + m 2 



1 + 



In 



1 + ^ + 1 



+ 



4C7f(1-0 



7 + 9£ + 2? 

s 2 — 2m 2 
s 2 + m 2 



2m 2 (l-i) 



1 + 



6-£ + 



- 1 
2m 2 (l 



1 , m 
+ - In — 

2 fj,' 







In 



4s 2 
H 2 



(4.1) 



Ca 9 + 7£ + 2£ 2 



6m 2 (l -£) 



s 2 + m 2 



(6 - i6e + e 2 



2m 2 



2m 4 (1 -£) 



(s 2 — m 2 ) 2 



[9 - 4£ + 



2m 2 (l -0 



' s 2 + m 2 m 2 ( s 2 
m = 1 =- m lH ^ 

We find that the nonperturbative T5 is several orders of magnitude larger than the one-loop 
form, and there is no sign of the lattice data approaching the perturbative form even for 
the most ultraviolet points we can trust, around 5 GeV. We take this as an indication that 
very strong nonperturbative effects affect this form factor. It is also worth noting that the 
one-loop contribution to both T5 and X[ at the symmetric point are an order of magnitude 
smaller than the one-loop contributions to form factors at the asymmetric point. 

In order to get a dimensionless measure of the strength of this component relative to 
the tree-level vertex, we have scaled T5 with the gluon momentum q. We show this together 
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Q(q) = 2K(p) (GeV) 

Figure 11: The dimcnsionless form factors A^ and qr§ at the symmetric point, as a function of q. 
These quantities gives a measure of the relative strength of the two components of the vertex. 

with A' x in figure |ll|. As we can see, between 1 and 2 GeV, T5 contributes with about the 
same strength as A' l5 making it a very significant contribution that cannot be ignored. 

Although T5 has the same tensor structure as the (chromo-)magnetic moment, the 
relation between the two is not straightforward. In particular, since quarks are never on- 
shell, the Gordon decomposition which is used to define the magnetic moment in QED is 
not applicable, making the definition of the chromomagnetic moment ambiguous. This is 
an issue that deserves further investigation. 

5. Outlook 

We have computed the complete quark-gluon vertex at two kinematical points, finding 
substantial deviations from the abelian form — which cannot be described by a universal 
function multiplying the abelian form as in ||]. This, and the fact that we observe a in- 
dependent enhancement of Ai at the asymmetric point, where q and thereby also the ghost 
form factor G(q 2 ) is fixed, indicates that the ghost-quark scattering kernel entering into 
the Slavnov-Taylor identity ( |1.1| ) must contain nontrivial structure. 

The form factor T5, related to the chromomagnetic moment, has been estimated non- 
perturbatively for the first time, and found to be important. The work has been carried out 
on a relatively small lattice, using a fermion discretisation which has serious discretisation 
errors at large momenta. It will be important to repeat this study using larger lattices and 
a more well-behaved fermion discretisation. 

A natural extension of this work would be to map out the entire kinematical space in 
the three variables p 2 ,k 2 ,q 2 . This is numerically very demanding, but work is underway 
on a complete determination of X[(p 2 , q 2 , k 2 ) . 

Finally, it should be noted that the lattice Landau gauge restricts us to computing 
only the transverse-projected vertex away from q = 0; i.e., it is not possible to determine 
Xi,Ti,i = 1,...,4 separately; only the linear combinations A^. Although the vertex is 
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always contracted with the gluon propagator in all actual applications, and thus only the 
transverse-projected vertex plays any role in Landau gauge, it would be of interest to 
determine all these form factors by computing the vertex in a general covariant gauge - 
which would also give a handle on the important issue of gauge dependence. 
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A. Tree- level expressions 



The tree- level lattice expressions are given in terms of the lattice momentum variables, 



2 / 

Qu(p) = -sm(p fl a/2) = — Jl- cos(p„a) ; (A.2) 
a a v 

K^p) = -K^p) = i- sin(2p M a) ; (A.3) 
C M (p) = cos(p M a) . (A.4) 

The tree-level vertex is @ 

(P, q) = CnS® (p)" 1 ^ (p)A$ (p, q)S® (p + g)S<°> (p + q)- 1 , (A.5) 

and 

(Sf)- 1 ^ = ^A V + B V ) , (A.6) 

where we have written 



c m = 1 + b q m ; (A.7) 

A v (p) = 2c' g D(p); (A.8) 

B V (p) = (c m ~ 2c' q M{p))D(p) ; (A.9) 

Di{p) = [A 2 v (p)K 2 (p) + Blip)] /D{p) ; (A.IO) 

D{p) = K\p) + M 2 (p); (A.ll) 

Mip) = m + ^Q 2 (p). (A.12) 



At q = we have 

0,0) = ^(t jrA v (p)+B v (p)) (7 M C M (p)-^(p)) (i ^A v (p)+S v (p)) . (A.13) 
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This expands to 



(i KA V + B V ) (7^ - (i KA V + V 

=(A 2 V K 2 + B V )^C^ - 2A 2 v C tx KK^ + 2iA v B v K^C fJ 
- iB v K^ + 2A v By KK^ + iA v K 2 K^ 

=(A 2 V K 2 + B V )C^ + 2{A V B V - A 2 V C^) KK^ 
+ i(2A v B v C^ + A 2 V K 2 - B V )K^ . 

The tree- level form factors Ajj j A3 can be read off directly: 



(A.14) 



x(0) , 1(0) r 

A 3 + a 3 o M 



CmC, 



,D 2 



q D 2 i 



c m - 2c' q M - 2c' q C^ 



2D 2 



" 1 AyK 2 — By + 2AyByCf, 



(A.15) 
(A.16) 



The lattice, tree- level corrected equivalents of ( 2.12 ) and (2^9), which we use to obtain A2 
and A3, are thus 



a 2 Cp 2 ,o, p 2 



— y 



AK 2 {p) 



1 



ImTr 7/ ,A M (p,0) + Ai(/,0,p 2 



Hp 2 Ap 2 



— y 



2K 2 (p) 



1 

4^o 



(A.17) 



^(p)-^ReTr A M (p,0) - 2K 2 (p)(xf ] (p) + Af(p)C„(p) 



(A.18) 

The tree-level vertex at the symmetric point is given by eq. (B.21) of 0. We use the 
following decomposition into independently transverse tensors, 



(°), 



fjo 



4 \K-QC^- KK, 



?(°)/ 



(A.19) 



~,(0) 



ir 5 



[c„ y (T^uQu + K^J2 °v\Q»K x /{Q ■ K) 



where Q = Q(q),K = K(q),C = C(q/2); X^\tP and f 3 w are given by (B.25)-(B.27) of 
H, and 



.(0) 



,(0) 



-,(0) 
^5 



c D 2 (n) 
j^A v {p)By(p) = 2c m c> q (c m - 2c' q M(p))^ ; 



c sw c m Djp) 
2 Diip) ' 
A 2 yjp)iKjp) ■ Kj P )) 
Djip) 



4c swCmC ' g 2 (K(p).K(p))^|| 



(A.20) 
(A.21) 
(A.22) 
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Since the continuum A' 1 (7 /J 



becomes two independent tensors on the lattice, 



\l I I 2 _2 I 2 

A i I T ,q 'T 



7m 



(Ai-Q 2 r 3 ) 



Q 2 



~h{K-QC ulu - pK^) , (A.23) 



we cannot simply factor out the tree-level behaviour with a simple multiplicative correction. 
Instead we apply a 'hybrid' scheme where the dominant term, multiplying (7^— QQ^/Q 2 ), 
is corrected multiplicatively, after first subtracting off the remaining part, 



{K-QC^- ^K u ) - K-QC^ 



JQ 2 



4 0) {K 2 -K,Q^). (A.24) 



It turns out that this term is completely negligible, but it has still been included in the 
correction. Thus, the lattice, tree-level corrected equivalent of ( |2.1l| ) which we use to 
compute X[, is 



2 q *\ l_-4 ImTr ^(-9Ag)-ff(^-^Q^) 



A 



(0) 



,(0), 



(A.25) 

For T5 we employ an additive correction scheme, and thus the lattice equivalent of ( 2.1C| ) 
is 



, q 2 2 1 2 



1 

'4<? 



Q M (g)^RelV^Af(-g/2,g) 



Ql{4° ] + C„C„f 5 (0) + C » ? ? [1 - {Cu - C^Ql/Q -K) } 



(A.26) 
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